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Abstract. We present results derived from a high-resolution cosmolog- 
ical iV-body simulation in which the equations of motion have been 
changed to account for Modified Newtonian Dynamics (MOND). It is 
shown that a low-fio MONDian model with an appropriate choice for 
the normalisation as of the primordial density fluctuations can lead to 
similar clustering properties at redshift z = as the commonly accepted 
(standard) ACDM model. However, such a model shows no significant 
. structures at high redshift with only very few objects present beyond 

' z > 3. For the current implementation of MOND density profiles of 

, gravitationally bound objects at z — can though still be fitted by the 

f*"^ ■ universal NFW profile. 

O 1 Introduction 

a: 

| Although the currently favoured ACDM model has proven to be remarkably suc- 
cessful on large scales (cf. Spergel et al. 2003), recent high- resolution TV-body 
simulations seem to be in contradiction with observation on sub-galactic scales: 
the CDM " crisis" is far from being over. Suggested solutions to this include 
H I the introduction of self-interactions into collisionless A-body simulations (e.g. 
Spergel & Steinhardt 2000), replacing cold dark matter with warm dark matter 
(e.g. Knebe et al. 2002) or non-standard modifications to an otherwise unper- 
turbed CDM power spectrum (e.g. bumpy power spectra, Little, Knebe & Islam 
2003). Some of the problems, as for instance the overabundance of satellites, can 
be resolved with such modifications but none of the proposed solutions have been 
able to rectify all shortcomings of ACDM simultaneously. Therefore, alternative 
solutions are unquestionably worthy of exploration, one of which is to abandon 
dark matter completely and to adopt the equations of MOdified Newtonian Dy- 
namics (MOND; Milgrom 1983; Bekenstein & Milgrom 1984). 



Astrophysikalisches Institut Potsdam, An der Sternwarte 16, 14482 Potsdam, Germany 



© EDP Sciences 2008 
DOI: (will be inserted later) 



Mass Profiles and Shapes of Cosmological Structures 



Table 1. Model parameters. In all cases a value for the Hubble parameter of h — 0.7 
was employed. 

label Op 06 A gf=° o-g orm g [cm/s 2 ] 

ACDM CL30 U\04 OJ (188 088 = 

OCBM 0.04 0.04 0.0 0.88 0.88 

OCBMond 0.04 0.04 0.0 0.92 0.40 1.2 xl0~ 8 



2 The Simulations 

We adapted our cosmological iV-body code AMIGA 1 (Knebe, Green & Binney 2001) 
to account for the effects of MOND in the following way. In an TV-body code one 
usually integrates the (comoving) equations of motion 

x = -? , p = (2.1) 
a z a 

which are completed by Poisson's equation 

V, • F poc (x) = -A.$(f) = -AnG{p{x) - p) . (2.2) 

In these equations x — r/a is the comoving coordinate, p the canonical mo- 
mentum, Va;- the divergence operator (A^ the Nabla operator) with respect to 
x and F pec (x) = — V < E>(x) the peculiar acceleration field in comoving coordinates. 
We now need to modify these (comoving) equations to account for MOND. 

Sanders & McGaugh (2002) showed that the relation between the Newtonian 
<7 prop and MONDian g^j° p acceleration (in proper coordinates) can be written as 




where we already faciliated Milgrom's interpolation function p(x) — x(l + x 2 )~ 1 / 2 
(Milgrom 1983) and go is the fundamental acceleration of the MOND theory. 

If we further assume that MOND only affects peculiar acceleration (in proper 
coordinates), i.e. g^,l° p = F pcc /a 2 , the recipe for adding the MOND formalism 
to an TV-body code reads as follows: (1) solve Eq. 1|2.2|) using AMIGA which gives 
the comoving F pec , (2) calculate the peculiar acceleration in proper coordinates 
gP™P = F pcc /a 2 , (3) use Eq. E3 to calculate f™* ec from gj^, (4) transfer g p ™* ec 

back to Fm, P cc — a2 5M°pec anc ^ (5) use -^M.pec rather than F pcc for the equations 
of motion l|2.1|l . For a more elaborate discussion of the assumptions upon which 
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this scheme is based and a more detailed derivation of the formulae presented here 
we refer the reader to Knebe & Gibson (2004). 

Our suite of simulations now consists of a standard ACDM model, an open, 
low-fio model with the same as as ACDM (OCBM), and an open, low-fio model 
with MOND and adjusted as (OCBMond), and their physical parameters are sum- 
marized in Tabled We simulated 128 3 particles in a box of side length 32ft -1 Mpc 
from redshift z = 50 to z = 0. Gravitationally bound objects were identified us- 
ing the AMIGA'S native halo finder AHF described in great detail in Gill, Knebe & 
Gibson (2004). 
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Fig. 1. Density field of the ACDM (left), OCBMond (middle), and OCBM (right) model 
at redshift z — 0. 



3 The Results 

The first intriguing result can be viewed in Fig. ^ where we show a projection of 
the whole simulation with each individual particle color-coded according to the 
local density at redshift z = 0. This figure demonstrates that the MOND model 
exhibits fairly similar features in terms of the locations of high density peaks, 
filaments and the large-scale structure when being compared to both the standard 
ACDM model and the OCBM simulation. One should bear in mind though that 
the OCBMond simulation was started with a much lower <7g orm normalisation than 
the other two run indicating a faster growth of structures in MOND universes. 

This result is supported by Fig. [2] (left panel) where the abundance evolution 
of gravitationally bound objects more massive than M > lO 11 ^ -1 M Q is shown. 
Moreover, this figure also poses a serious problem for cosmological MONDian 
structure formation. Due to the increased formation rate of objects we had to 
lower the amplitude of the primordial density perturbations. This in turn leads 
to a strong deficiency of gravitionally bound structures at redshifts z > 3, quite 
in contrast to observations where we find galaxies out to redshifts of z ~ 6 (Sh- 
ioya et al. 2005). 
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Fig. [3 (right panel) also shows the spherically averaged density profile of the 
most massive halo in all three models along with fits (thin solid lines) to a Navarro, 
Frenk & White (1997) profile p(r) oc ((r/r s )(l + r/r s ) 2 )~ 1 . We observe that even 
for the OCBMond model the data is equally well described by the functional form 
of a NFW profile (at least out to the virial radius indicated by the respective 
vertical lines). However, the central density of that halo in the OCBMond model 
is lower than in ACDM and especially in OCBM. A quantitative analysis further 
shows that the concentration of the OCBMond halo is about a factor of four 
smaller than in the ACDM model. 



Fig. 2. Redshift evolution of the abundance of objects with mass M > lO 11 /! -1 M 
(left) and the density profile for the most massive halo in all three models (right). 



4 The Conclusions 

Even though it is possible to match a cosmological simulation including the ef- 
fects of MOND to the standard ACDM structure formation scenario at redshift 
z = there are serious deviations at higher redshifts. We conclude that the most 
distinctive feature of a MONDian universe is the late epoch of galaxy formation. 
However, the density profiles of gravitationally bound objects still follow the uni- 
versal NFW shape. 
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